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Abstract 

For the quantum kicked top we study numerically the distribution of Hilbert- 
space vectors evolving in the presence of a small random perturbation. For an 
initial coherent state centered in a chaotic region of the classical dynamics, the 
evolved perturbed vectors are distributed essentially like random vectors in 
Hilbert space. In contrast, for an initial coherent state centered near an elliptic 
(regular) fixed point of the classical dynamics, the evolved perturbed vectors 
remain close together, explore only few dimensions of Hilbert space, and do 
not explore them randomly. These results support and extend results of earlier 
studies, thereby providing additional support for a new characterization of 
quantum chaos that uses concepts from information theory. 

I. INTRODUCTION 

In a series of previous papers two of the authors introduced a characterization of 

Hamiltonian chaos which is directly applicable to quantum as well as to classical systems. 
This characterization, formulated in the framework of statistical mechanics, is based on the 
following question: How much information is needed to predict an evolved system state 
in the presence of random perturbations? General arguments p| and an investigation of 
the symbolic dynamics of the baker's map provide strong evidence that chaotic classi- 
cal Hamiltonian systems show what we call hypersensitivity to perturbation — i.e., a rapid 
increase with the number of time steps of the information needed to describe the perturbed 
time evolution of a system state, the information attaining values exponentially larger than 
the increase of ordinary entropy that results from averaging over the perturbation. 

Hypersensitivity to perturbation explains quantitatively how information about the state 
of a system is lost through interaction with an incompletely known environment and, there- 
fore, is important for understanding why entropy necessarily increases in systems that are 
not perfectly isolated. The connection of this work on chaos with statistical physics is devel- 
oped in Section |V|. In addition to its straightforward motivation in statistical mechanics. 



the concept of hypersensitivity to perturbation may provide a more physical way to charac- 
terize quantum chaos f^. Numerical simulations show that the quantum baker's map |^ 
displays hypersensitivity to perturbation. 

In this paper, we analyze how hypersensitivity to perturbation arises in the quantum 
kicked top a system whose classical dynamics has both chaotic and regular regions. To 
shed light on the reason for the rapid increase of information associated with the property of 
hypersensitivity to perturbation, we perform a detailed numerical analysis of how the vectors 
arising from different perturbation histories (realizations of the random perturbation) are 
distributed in Hilbert space. For an initial coherent state centered in a chaotic region of 
the classical dynamics, the evolved perturbed vectors are distributed essentially like random 
vectors in Hilbert space. In contrast, for an initial coherent state centered near an elliptic 
(regular) fixed point of the classical dynamics, the evolved perturbed vectors remain close 
together, explore only few dimensions of Hilbert space, and do not explore those dimensions 
randomly. 

Quantum systems show no "sensitivity to initial conditions," due to unitarity, but they 
show what one might call sensitivity to parameters in the Hamiltonian, as has been demon- 
strated for the kicked top by Peres 0. Peres compares the time evolution of the same initial 
Hilbert-space vector for two slightly different values of the twist parameter in the kicked-top 
Hamiltonian (see Section He finds that, after a fixed number of time steps, the two 
evolved vectors are far apart if the initial vector is a coherent state centered in a chaotic 
region of the classical dynamics, but the two evolved vectors stay close together if the initial 
coherent state is centered near an elliptic fixed point of the classical dynamics. 

Our approach to quantum chaos could be viewed as a generalization of Peres's work: 
while Peres studies time evolution due to an incompletely known Hamiltonian, we analyze 
the distribution of vectors arising from time evolution under a stochastic Hamiltonian. There 
is, however, a fundamental difference in philosophy between the two approaches, which can 
be understood fully only in the context of statistical mechanics. This difference becomes 
apparent in Section ^I[ 

The quantum kicked top, its classical limit, and the concept of a coherent state are re- 
viewed in Section Section |T| defines the perturbations used in the numerical simulations. 



In Section ^ we explain how the distribution of vectors in Hilbert space is connected to 
questions of information and entropy and how the numerical data are compiled into figures. 
Section contains the numerical results of this paper. Finally, in Section we discuss 
the implications of our results for the foundations of statistical physics. 



II. THE KICKED TOP 

The quantum model of the kicked top ||^J^ describes a spin- J particle — i.e., an angular 
momentum vector hJ = Ti{Jxi Jy, Jz)-, where [Jj, Jj\ = ieijkJk — whose dynamics in (2J + 1)- 
dimensional Hilbert space is governed by the Hamiltonian 

+ 00 

H{t) = {hp/T)J, + {nk/2J)Jl ^(^ - nT) . (2.1) 

n=— oo 

The free precession of the spin around the z axis (first term in the Hamiltonian) is interrupted 
periodically by sudden kicks or twists at times nT with twist parameter k (second term in 
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the Hamiltonian) . The angle of free precession between kicks is given by p. In this paper 
we always use p = 7t/2. 

We look at the time evolution of an initial Hilbert-space vector at discrete times 
nT. After n time steps, the evolved vector is given by 

m = mi^o) , (2.2) 

where Uk is the unitary Floquet operator: 

Uk = exp [-i{k/2J)J^] exp {-i7rJj2) . (2.3) 

The classical Poincare map corresponding to the quantum map is obtained by introducing 
the unit vector Co = {X, Y, Z) = 3/J and performing the limit J oo. One obtains [0 

X' = -r , 

Y' = XcoskY + ZsinkY , (2.4) 
Z' = ZcoskY - XsinkY . 



The map ( |2.4| ) is an area-preserving map of the unit sphere, i.e., an area-preserving map on 
the configuration space of a classical spin with fixed magnitude. Depending on the value of 
the twist parameter k, this map has regions of chaotic behavior interspersed with regular 
regions associated with elliptic cyclic points. In this paper, we are interested in two cyclic 
points of the map for k = 3. One is an elliptic fixed point of period 1 located at 

Z = cos e = 0.455719 , = 37r/4 , (2.5) 

where we have used spherical coordinates 

9 = arccos Z , 

V9 = arctan F/X . (2.6) 

The elliptic fixed point (|2.5| ) is surrounded by an oval-shaped regular region, extending 
about 0.3 radians in the (/^-direction and about 0.5 radians in the 6'-direction. The other 
cyclic point of interest to us is a hyperbolic fixed point of period 4, which has a positive 
Lyapunov exponent. It is located in the middle of a chaotic region at ||^ 

Z = cos9 = 0.615950 , = 7r/4 . (2.7) 



We choose the two fixed points (|2.5| ) and (|2.7| ) instead of the extreme elliptic {(p = rni/A 



Z = 0) and hyperbolic {Z = ±1) points described in Ref. |0, because the latter suffer 
accidental invariance with respect to one or both of the perturbation operators considered 
in Section |T|. 

The model has a conserved parity, which for half-integer J takes the form 

5 = — i exp(— zvr J^) (2.8) 

and which permits factorization of the matrix representation of the operator U into two 
blocks. Starting from a state with definite parity, the whole dynamical evolution occurs in 
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the invariant Hilbert subspace with the given parity. For half-integer J, the dimension of 
the even-parity subspace (eigenspace of 5* with eigenvalue 1) is J -|- |. In this paper, we 
work with J = 511.5 in the 5 12- dimensional even-parity subspace. We consider only the 
projection of the initial vector in the even subspace. Numerical evidence and symmetry 
considerations suggest that no additional insight is gained by including the odd-parity 
subspace. In any case, the restricted model can be regarded as a quantum map in its own 
right, which can be investigated independently of the behavior of the complete kicked-top 
model. 

We want to choose initial vectors for the quantum evolution that correspond as closely as 
possible to the classical directions ( |2.5| ) and (|2.7|). For this purpose, coherent states [p!0|-pr^ 
are appropriate. The coherent state \9,ip) is defined by the relation 

n- J\e,ip) = J\9,ifi) , (2.9) 

where n is the unit vector pointing in the direction given by 6 and (p. All coherent states can 
be generated by an appropriate rotation of the state \J,M = J) = \9 = 7r/2, ip = 0), where 
I J, M) (M = — J, . . . , J) is the common eigenstate of and with eigenvalues J(J + 1) 
and M, respectively. In calculations, it is convenient to use the explicit representation 

2J 



\e,ip) = J2^PA0,n)e''^^\J,n-J) , (2.10) 

n=0 



where 



PA9,n) = (^l±^p^y'"" l^ l-.rccos9 j ^^^^^^ 

In the following we need a metric on Hilbert space. The distance between two normalized 
vectors lipi) and \^|J2) is defined as the Hilbert-space angle 

s(|V^i),|V^2)) = cos-i(KV^i|^2)|) = (2.12) 

between the two vectors [13]. Consider two coherent states \6,(p) and \6',ip'). In terms of 
the angle a between the directions (6*, ip) and {9', ip'), the distance between the two coherent 
states is given by [|14 



cos[3(|^^,^),|^^',^'))] = I e',^')\ = [cos(«/2)]^-^ exp(-JaV4) , (2.13) 

where the approximation is valid for large J. Two coherent states can therefore be regarded 
as roug hly orthogonal if a > 2J-^/^ [g. The size of the coherent state \6, ip) is conveniently 
defined in terms of the Q function 

QeA(^',ip') = \{9',ip' I e,ip)\' = [cos(a/2)]^-^ ^ Q(«) . (2.14) 

Since Q(2 J^^/^) ^ e^'^Q{0), the Q function of the coherent state \9, ip) is very small outside 
a region of radius 2J^^/^ centered at the direction {9,ip). For the value J = 511.5 used 
in this paper, one finds a radius of 2J~^/^ ~ 0.09 radians, less than the size of the regular 



region around the elliptic fixed point ( p.5|) 
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III. PERTURBED EVOLUTION 



Our goal is to quantify how much information is required to track the state of a system 
when, instead of being wholly isolated, it is perturbed by interaction with its environment. 
In classical physics, interactions with an incompletely known environment can be described 
by a stochastic Hamiltonian, each realization of which corresponds to a particular initial 
condition for the environment. The situation is more complicated in quantum mechanics. 
Due to the possible entanglement of system and environment, there is no way to associate a 
unique perturbation history with a given initial condition of the environment. An upcoming 
publication discusses how to study hypersensitivity to perturbation in a realistic model 
of a quantum system interacting with an environment. For the present paper we restrict 
ourselves to the special case of quantum time evolution under a stochastic Hamiltonian. 

The problem is simplified further by considering only two possible unitary time evolutions 
at each step. These two time evolution operators we denote by If^ and U~ . A perturbed 
time step consists in applying either f/~*" or U~ with equal probability. After n time steps, 
the number of different perturbation sequences is 2", each sequence having probability 2~". 

We use two different perturbations: (i) the twist perturbation, defined by choosing the 
twist parameter k at random from step to step, the two possible Floquet operators being 
given by §] 

u+ = Uk, u- = Uk+e; (3.1) 

(ii) the turn perturbation, defined by rotating the spin by a small angle e around the z axis 
after each unperturbed step Uk, the two possible Floquet operators being given by 

f/+ = exp{-ieJ,) Uk , = exp{+ieJ,) Uk ■ (3.2) 



Notice that the time-evolution operators (O) and (^) commute with parity ( f^.8| ) and 



hence do not couple odd- and even-parity subspaces. 

For perturbation strengths we use e = 0.03 and e = 0.003 for the twist perturbation 
and e = 0.003 for the turn perturbation. For a twist parameter + e the zenithal location 
of the fixed points and (^) changes slightly to Z = 0.443579 and Z = 0.619848 for 
e = 0.03, and to Z = 0.454497 and Z = 0.616341 for e = 0.003. This corresponds to changes 
in the zenithal angle 6 of ~ 0.014 radians and AO ~ 0.005 radians for e = 0.03, and to 

~ 0.0014 radians and A^ ~ 0.0005 radians for e = 0.003. All these angles, as well as 
the angle e = 0.003 we use for the turn perturbation, are very small compared to the size of 
the elliptic region around the fixed point ( p.5| ) and are also much smaller than the size of a 
coherent state. 



IV. DISTRIBUTION OF VECTORS AND INFORMATION 

The 2"- different perturbation sequences obtained by applying every possible sequence of 
tj~ and for n time steps lead to a list of 2" vectors, each having probability 2~". In 
this section, we explain how the distribution of these 2" vectors in Hilbert space is related 
to information and entropy. 
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Let us start with a slightly more general situation. Imagine we are given a list of N vec- 
tors in D-dimensional Hilbert space, {ipi), ■ ■ ■ , \iPn), with probabilities Pi, ■ ■ ■ ,Pn- Together 
with our knowledge of the system Hamiltonian and boundary conditions, the list of vectors 
with their probabilities constitutes our background information. We ask for the average in- 
formation needed to specify a single one of these vectors, given the background information. 
The information to specify a particular vector can be quantified either via conditional algo- 
rithmic information or by the length of a codeword in some coding scheme |]TB|]. In both 



cases, it is a consequence of the variable-length coding theorem [|T6| that the information 



averaged over all vectors, or average information, is bounded below by 

N 

A/= -^paOg2P^ . (4.1) 

1=1 

(Throughout this paper, information and entropy are measured in bits.) There exist coding 



schemes — an example is Huffman coding ||T^ — where AJ + 1 is an upper bound for the 



average information or codeword length. It can be shown |T^ that there exists a universal 
computer for which AI -|- 1 is an upper bound for the average algorithmic information as 
well. Therefore, we call AI the average information, in the sense that the actual average 
information is within one bit of AI if an efficient coding scheme is used. In the case of 2"' 
equiprobable vectors the average information is AJ = n. 

Suppose some of the vectors are very close together in Hilbert space, so that they 
form a small group. If one is interested in lowering the amount of information AI, one may 
choose to provide just enough information to specify that the actual vector is located in 
that group, the price being that the entropy of the group is generally bigger than zero. To 
be more specific, we introduce a coarse graining on Hilbert space defined by a resolution 
angle 0. Vectors less than an angle apart are grouped together. More precisely, groups 
are formed in the following way. Starting with the first vector in the list, I'^'i), the first 
group is formed of and of all vectors in the list that are within an angle of The 
same procedure is then repeated with the remaining vectors to form the second group, then 
the third group, continuing until no ungrouped vectors are left. This grouping of vectors 
corresponds to a partial averaging over the perturbations. To describe a vector at resolution 
level amounts to averaging over those details of the perturbation that do not change the 
final vector by more than an angle 0. 

For a given resolution 0, there are A^(0) groups. We denote by Nj the number of vectors 
in the jth group 

{Ef=t^ Nj = N). The Nj vectors in the jth group and their probabilities are 
denoted by \ipi), . . . , lipj^.) and pi, ■ ■ ■ iP'njI respectively. Knowing that the system state is in 
the jth group, but not knowing which state in the jth group is the actual state, corresponds 
to describing the system by the density operator 

h = <u'T.pi\^r)^i\^ (4-2) 

where 

Nj 

q^ = Y.Pi (4-3) 

i=l 
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is the probability of the j'th group given only the background information. For = 7r/2, there 
is only one group, whose density operator, denoted by p{n/2), corresponds to a complete 
average over the perturbations. 

The average information needed to specify which group a vector is in — i.e., the average 
information needed to specify the system state at resolution level — is given by 

N{<l>) 

A/(0) = -^g,log2g, . (4.4) 

The entropy of the jth group is given by the von Neumann entropy 

Ai7^ = -Tr[/3,log2p,] , (4.5) 
and the average entropy, called trade-off entropy in the following, is 

AH{(f)) = J2 lAH^ ■ (4-6) 

i=i 

If one chooses to describe the set of vectors not exactly, but only at resolution level 0, 
the average information needed to specify the system state decreases. There is, however, a 
price: with increasing resolution angle, the uncertainty about the system state increases on 
the average, to a degree quantified by the trade-off entropy. 

As a further characterization of the way the vectors are distributed in Hilbert space, we 
want to define a quantity that indicates how many dimensions of Hilbert space are explored 
by the vectors in a group. One such quantity would be the dimension of the subspace 
spanned by the vectors, which is equal to the number of nonzero eigenvalues of pj. In 
practice, however, this is not a very useful measure because it cannot discriminate between 
the case in which all dimensions are occupied with equal weight (all eigenvalues of pj roughly 
equal) and the case in which most vectors are concentrated in a low- dimensional subspace 
(all eigenvalues of pj nonzero, but of strongly varying magnitude). 

A possible measure of the number of explored dimensions, which takes into account the 
small weight of dimensions corresponding to relatively small eigenvalues, is the exponential of 
the entropy, 2^^^ . This quantity is bounded above by Dj if the vectors are confined to a Dj- 
dimensional subspace and gets smaller if the dimensions are occupied with different weights. 
For example, if two eigenvalues of pj are close to | and all the others are close to zero, then 
2^H^ ~ 2, indicating that the vectors are essentially confined to a two-dimensional subspace. 
Unfortunately, for small resolution angles 0, Aif-' is necessarily small just because all the 
vectors in the group lie along roughly the same direction in Hilbert space; this is true even 
if the orthogonal components of the vectors are evenly distributed over all the orthogonal 
directions in Hilbert space. For example, the density operator describing a uniform distribu- 
tion of vectors within a sphere of radius ^ 7r/2 has one dominating eigenvalue close to 1 
and D — 1 eigenvalues that are all equal and close to zero (see Appendix P) . Clearly, in this 
case 2^^^ is not an adequate measure of the number of dimensions explored. On the other 
hand, if one could disregard the largest eigenvalue in this example, then the exponential of 
the entropy would still be a useful measure of the number of explored dimensions. 

We therefore introduce the spread AH2 as the entropy calculated with the largest eigen- 
value of Pj omitted. The spread is defined as 
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where A{ > A2 > • • • > are the nonzero eigenvalues of the density operator pj. The 
average spread is 

AH2i<P) = E lAHi ■ (4.8) 

i=i 

For a given resolution angle 0, the entropy AH^ is bounded above by the entropy -f/^D,max(0) 
of Eq. ( P10|) , which for small has the value i?D,max(0) — 0^ ^og2[e{D — l)/0^]. The spread 
AH2, on the other hand, can attain its maximum value AH2 = log2(-D — 1) for arbitrary 
resolution angles 0. Indeed, the entropy and the spread are related by 

AH^ = -Xi log2 Ai - (1 - Xi) log2(l - Aj) + (1 - Xi)AHi , (4.9) 

which indicates how, when there is one dominating eigenvalue, a large spread does not lead 
to a large entropy. 

By giving different weight to dimensions corresponding to different eigenvalues of pj, the 
quantity [2^^2] turns out to be a good indicator of the number of Hilbert-space dimensions 
explored by the vectors in a group, independent of the size of the region occupied by the 
group, ([x] denotes the smallest integer greater than or equal to x.) In our analysis of 
the numerical results, we identify the number of dimensions explored by the total set of 
vectors with the integer Ud = ^2^-^'^^'^^'^'']. 

By determining the information A/(0), the trade-off entropy Aif (0), and the average 
spread Ai72(0) as functions of the resolution angle 0, a rather detailed picture emerges of how 
the vectors are distributed in Hilbert space. The information summarizes the distribution 
of group sizes at the given resolution. The trade-off entropy and the average spread indicate 
how the vectors are distributed inside the groups. 

It is easy to see that information and trade-off entropy obey the inequalities 

A/(0) > AJ(0) > AJ(7r/2) = (4.10) 

and 

= AH{0) < Ai7(0) < AH{tt/2) . (4.11) 

The first inequality in ( [4.10| ) follows from the fact that any group at resolution is the 
union of groups at resolution = 0; in words, the average information needed to specify a 
group at resolution = is equal to the average information needed to specify a group at 
resolution plus the average information needed specify = groups within the groups at 
resolution 0. The last inequality in Eq. ( [4.11| ) is a consequence of the concavity of the von 



Neumann entropy |T9|,^. A general theorem about average density operators |T9|,^ shows 
that, for all 0, 

A/(0) + AH{(f)) > AH{n/2) . (4.12) 
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In general, A/(0) is a decreasing function of 0, whereas AH{(f)) is increasing. This mono- 
tonicity can sometimes be violated, however, because of discontinuous changes in the group- 
ing of vectors. 

As a still further characterization of our list of vectors, we calculate the distribution g{(j)) 
of Hilbert-space angles (p = s{\ip), \ip')) = cos^^^Kiplip')]) between all pairs of vectors \ip) 
and For vectors distributed randomly in Z)-dimensional Hilbert space, the distribution 
function g{(j)) is computed in Appendix ^ 

g{<P) = ^^^MI^ = 2{D - I)(sin0)2^-=^cos0 (4.13) 



[Eq. ( [A24D ]. Here Vd(0) = (sin 0)2(^"^)Vd [Eq. ( |AT8| )] is the volume contained within a 



sphere of radius in D-dimensional Hilbert space, and Vd = vr^^^/(D — 1)! [Eq. ([A22D ] is 
the total volume of Hilbert space. The maximum of g{(f)) is located at = arctan(A/2i5 — 3); 
for large-dimensional Hilbert spaces, g{(j)) is very strongly peaked near the maximum, which 
is located at ~ 71/2 — l/y/2D — 3, very near 71/2 (see Fig. |1]). 



V. NUMERICAL RESULTS 

In this section, we describe the numerical results, which are shown in the figures. In all 
numerical examples, we use spin J = 511.5 and unperturbed twist parameter k = 3. The 
calculations are done in the 512-dimensional even-parity subspace (eigenspace of the par- 
ity ( |2.8| ) with eigenvalue 1), i.e., effectively in a 512-dimensional Hilbert space. Throughout 
this section, we use only two different initial states. The first one, the coherent state \6,(p) 
with 6 and cp given by Eq. ( |2.5| ), is centered in a regular region of the classical dynamics; 
we refer to it as the regular initial state. The second one, referred to as the chaotic initial 
state, is the coherent state \6, if) with 9 and given by Eq. (|2.7|) ; the chaotic initial state is 
centered in a chaotic region of the classical dynamics. Subsections |V A| and V B describe re- 



sults for the twist perturbation with two different perturbations strengths e. Subsection [VC 
contains results for the turn perturbation. 



A. Twist perturbation: e = 0.03 

Applying all possible n-step perturbation sequences, i.e., all possible sequences of n 
Floquet operators f/" and to an initial state iV'o) generates a set of 2" equally likely 
states, as considered in Section In Fig. |], the quantities defined in Section are 
computed for the twist perturbation ( p.l|) with perturbation strength e = 0.03. Figure |^ 
shows results after 8 and 12 steps for both the chaotic and the regular initial conditions. 

In the chaotic case in Fig. 0(a), where the total number of vectors is 2^ = 256 (8 steps), 
the distribution of Hilbert-space angles, gi^cj)), is concentrated at large angles, i.e., most pairs 
of vectors are far apart from each other. The small peak of g{(f)) at ~ 7r/16 corresponds to 
128 pairs of vectors, the two vectors in each pair being generated by perturbation sequences 
that differ only at the first step. The somewhat larger peak of g{(j)) at ~ 37r/16 similarly 
indicates the existence of 64 quartets of vectors, generated by perturbation sequences dif- 
fering only in the first two steps. The information AJ needed to track a perturbed vector 
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at resolution level is 8 bits at small angles where each group contains only one vector. At 
~ 7r/16 the information drops to 7 bits, and at ~ Svr/IG it drops to 6 bits, reflecting 
the grouping of the vectors in pairs and quartets, respectively. For larger resolution angles, 
the information stays constant before dropping rapidly to zero at angles ^ Svr/S. Just 
as the information begins to drop rapidly, there is a sudden drop to about 5 bits, reflect- 
ing a further, approximate grouping into 32 octets of vectors, generated by perturbation 
sequences that differ only in the first three steps. The final drop in the information coin- 
cides with the main peak in the angle distribution g{(j)) and with the rising of the trade-off 
entropy to its maximum value of AH ~ 4 bits. The number of explored dimensions is 

If the number of steps is increased to 12 [see Fig. |^(c)], the main features of Fig. ||(a) 
are preserved. The discontinuous drops in information — from 12 to 11 and from 11 to 10 
bits — due to the formation of pairs and quartets are obvious, but the corresponding peaks 
in g{(j)) are now almost invisible due to the larger scale produced by the larger total number 
of pairs of vectors. There is now little evidence of a further grouping into octets. Indeed, 
Figure ^(c) suggests that, apart from the organization into pairs and quartets, there is not 
much structure in the distribution of vectors for a chaotic initial state. The 1024 quartets 
seem to be rather uniformly distributed in a = \2^-^^^'^^'^^~\ = 65-dimensional Hilbert 
space. In order to check that the quartets are indeed more or less randomly distributed in 
Hilbert space. Fig. ^(c) should be compared to Fig. 0(b), where 1024 random vectors in a 
62-dimensional Hilbert space are shown. The random vectors are chosen at random from an 
ensemble distributed uniformly over Hilbert space The number of dimensions, 62, was 
chosen by trial and error so that the total entropy, Aif(7r/2), came out to be the same in 
Figures ^(c) and |l|(b). Although the angles between the random vectors are concentrated 
somewhat more towards larger angles, there is a striking similarity between these two figures. 

As further evidence of the nearly random character of the distribution in Fig. |^(c). Fig. |^ 
compares the eigenvalues of the density operators p(7r/2) corresponding to Fig. ^(c) and to 
the random vectors in Fig. |l|(b). The 62 largest eigenvalues in the chaotic case are almost 
identical to the 62 eigenvalues corresponding to random vectors in 62-dimensional Hilbert 
space. 

The distribution of perturbed vectors starting from the regular initial state is completely 
different from the chaotic case. Figure 0(b) shows the regular case after eight steps. The 
angle distribution g{(j)) is conspicuously non-random: it is concentrated at angles smaller 
than roughly 7r/4, and there is a regular structure of peaks and valleys. The information 
drops rapidly, with little plateaus corresponding to the valleys in the angle distribution. The 
number of explored dimensions is = 2, which agrees with results of Peres |Q that show 
that the quantum evolution in a regular region of the kicked top is essentially confined to a 
2-dimensional subspace. 

Figure 0(d) shows the regular case after 12 steps. The average information and the trade- 
off entropy show a behavior similar to the 8-step case. The plateaus in the information are 
washed out, corresponding to less pronounced minima in the distribution of angles. This 
appearance of formerly forbidden angles is expected as the number of vectors increases; it 
would occur even for a completely regular array of vectors. The eigenvalues of the density 
operator /)(7r/2) corresponding to the 4096 vectors in Fig. 0(d), shown in Fig. 0, confirm the 
restriction to a 2-dimensional subspace. The third-largest eigenvalue, measuring the relative 



10 



weight of the third explored dimension, is of order only 10~^. 

The results shown in Figures and |^ display a striking difference in the distribution of 
vectors in the chaotic and regular cases. In the chaotic case, the vectors, aside from the 
quartet structure, are distributed randomly in a subspace whose dimensionality increases 
with the number of steps. The information needed to track a perturbed vector after n steps 
is of the order of n bits, similar to the information needed to specify a vector out of a set 
of 2" random vectors. By contrast, in the regular case the vectors do not get far apart in 
Hilbert space, explore only few dimensions, and do not explore them randomly. 



B. Twist perturbation: e = 0.003 

Figure ^ shows the distribution of vectors arising from perturbed evolution with a very 
small twist perturbation of strength e = 0.003. Here, after 12 steps the vectors are not 
spread very far in Hilbert space. This is true even in the chaotic case, shown in Figure ^(a), 
where a typical angle between vectors is = 7r/8, the information decreases rapidly with the 
resolution angle, and only = 6 dimensions are explored. Even so, the chaotic case can be 
easily distinguished from the regular case, shown in Fig. §(b), where the perturbation has 
almost no effect on the time evolution of the vectors. 

To get a picture of the distribution of vectors for a larger number of steps. Fig. ^(c) 
shows 4096 vectors selected randomly from the 2^*^° vectors after 200 steps in the chaotic 
case, and Fig. ^(d) shows 1024 vectors selected randomly from the 2^*^° vectors after 200 
steps in the regular case. In the chaotic case, the 4096 vectors fill an = 373-dimensional 
subspace quasi-randomly, as can be checked by comparison with Fig. |T](a), where results for 
4096 random vectors in 512 dimensions are shown. In the regular case, shown in Fig. §(d), 
even after 200 steps not more than = 2 dimensions are explored. The vectors remain 
very close together, and the information drops rapidly with increasing resolution angle. The 
difference between the chaotic and regular cases is as striking as in the previous subsection. 

Although the data shown in Fig. ^(c) establish that the 4096 vectors selected from the 
available 2^°° vectors in the chaotic case fill a large-dimensional space quasi-randomly, they 
by no means establish that the distribution of all 2^°° vectors is similar to the distribution 
of 2^*^° random vectors. For example. Fig. ^(c) would look exactly the same whether the 
2200 vectors were randomly distributed or were organized into 2^"°""* randomly distributed 
groups, each consisting of 2"^ tightly bunched vectors, provided that the probability to select 
more than one vector from any group is negligible — i.e., provided that (4096)^/2^'^°""^ = 
2m-i76 ^ -^^ which is satisfied for m < 170. Indeed, in view of the data in the previous 
subsection, one expects the vectors to be organized on small angular scales into pairs and 
quartets and perhaps into somewhat larger groups that persist from the first few steps. 

To characterize the angle distribution completely would require the computation of the 
angles between all pairs among the 2^°° vectors, which would exhaust the storage and com- 
puting power of any computer now and in the foreseeable future. Our results are thus 
rigorous only up to 12 steps, where we are able to compute the angles between all pairs of 
vectors. Nonetheless, our results provide some support for the conjecture that the distribu- 
tion in the chaotic case is essentially random for large numbers of steps. In order to give 
an approximate picture of such a random distribution, we have developed approximations. 
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shown in Fig. ^ for the information and the trade-off entropy for 2^°° random vectors in a 
5 12- dimensional Hilbert space. 

These approximations are based on knowing A'^d jnax(0); the maximum number of disjoint 
spheres of radius that Z)-dimensional Hilbert space can accommodate. In Appendix |^ 
we compute the volume Vd(0) = {sm(f)y^^^^^VD contained within a sphere of radius cj) 
in Z)-dimensional Hilbert space [Eq. (|A18|) ] and the total volume Vd = ti^^^/ID — 1)! of 



D-dimensional Hilbert space [Eq. ( [A22[ )], from which it follows that 



A^D,max(0) = = (sin0)-2(^-^) . (5.1) 

It is worth emphasizing just how enormous Hilbert space is by noting that the number of 
spheres of radius (p = 0.1 rad that can be accommodated within a D-dimensional Hilbert 
space is A'^o.max — 0"^^^"^-* = 10^^^"^^; for the 512-dimensional Hilbert space considered in 
this paper, this is 10^°^^ spheres. 

Suppose A^ vectors, distributed randomly in D-dimensional Hilbert space, are grouped 
at resolution level (p. The number of groups at this resolution, N{(f)), cannot be larger than 
A^, the total number of vectors, nor larger than A^^d jnax(0); the maximum number of groups. 
For the average information, this entails 

A/(0) < log2 Ar(0) < AlD,m.M = min(log2 ND,m.M Aog^ n) . (5.2) 

One expects a large number of random vectors to fill Hilbert space almost uniformly. For 
angles for which A''^) max(0) ^ A^, there are close to ND^^g^^^cp) groups with roughly equal 
numbers of vectors in each group; therefore, for those angles 0, AI{(f)) ~ logg A?"/) max(0)- 
For angles for which A^D,max(0) ^ A^, there is just one vector in each group, whence 
AJ(0) ~ log2 A^ for those angles. This means that the upper bound AJ£) inax(0) is an 
excellent approximation to A/(0) everywhere except for a small region near the sharp bend 
located at the angle 06 determined by ND,uia.yi{,<t'b) = N. The upper bound AJ£)_max(0) is 
plotted in Fig. | for D = 512 and A^ = 2^00^ 

The trade-off entropy, on the other hand, cannot be larger than the maximum entropy 
of a group, i.e., 

AH{<f)) < HD,mM) - Hd{<P) , (5.3) 

where -ffi:),max(0) [Eq- ( P10|) ] is the maximum possible entropy for a density operator con- 
structed from vectors that lie within a sphere of radius in D-dimensional Hilbert space, 
and if£)(0) [Eq. (|B6|)] is the entropy of a uniform distribution of vectors within a sphere of 
radius 0. For large-dimensional Hilbert spaces, there is no appreciable difference between 
-f^D,max(0) and H£){(j)). The maximum entropy Hj:,^^i^{.(p) is plotted in Fig. |^ for D = 512. 

For angles for which A'^c),max(0) ^ A^, where the number of vectors per group is large 
enough — say, larger than D — the distribution of vectors within each group approximates 
a uniform distribution, and thus AH{(f)) is well approximated by H £>{(!)) and, hence, by 
-f^L>,max(0)- This mcaus that there is a region to the right of (pb in Fig. |^, where Hri^^s^^{(p) 
is not only an upper bound, but is also a good approximation to the trade-off entropy. For 
angles to the left of 06, the number of vectors in a typical group rapidly approaches 1, 
which means that the trade-off entropy is very close to zero in the region where the average 
information saturates. 
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C. Turn perturbation: e = 0.003 



Figure ^ displays results for the turn perturbation, showing the same range of behavior as 
in the preceding subsections. Figure ^(a) shows how the 2^^ vectors generated after 12 per- 
turbed steps in the chaotic case fill an ra^ = 46-dimensional Hilbert space randomly, except 
for a grouping into pairs, quartets, and perhaps octets, corresponding to the discontinuous 
drops in AJ(0). By contrast, in the regular case after 12 steps, displayed in Fig. ||(b), the 
vectors remain close together and fill just = 2 dimensions. Figures ^(c) and (d) show 2^° 
vectors chosen randomly out of the 2^° vectors generated after 30 perturbed steps. In the 
chaotic case, the distribution can barely be distinguished from the distribution of random 
vectors in Fig. |l](a). In the regular case, the vectors are spread a little further apart in 
comparison with Fig. |^(b), but they still fill only = 2 dimensions. 

Our results establish well the nearly random character of the distribution of vectors in 
the chaotic case. This is the main result of this paper, providing numerical evidence for 
hypersensitivity to perturbations in the quantum kicked top. It is more difficult (and less 
interesting) to give a general characterization of the distribution of vectors in the regular case. 
One reason for this is the finite size of the regular region on the classical unit sphere, which 
makes possible a sort of diffusion of a perturbed vector out of the regular region into the 
chaotic region. Figure |^ investigates this kind of behavior by showing the distribution of 2^° 
vectors randomly chosen after 100 [Fig. 0(a)] and 200 [Fig. ^h)] steps in the regular case. It 
is apparent that the vectors drift more and more apart and begin to explore more dimensions 
of Hilbert space, although even after 200 steps, the number of explored dimensions is still 
only rirf = 5. 

Figure |^(c) shows the eigenvalues of the density operators obtained by averaging over 
2^^ vectors randomly chosen after 30, 100, and 200 perturbed steps in the regular case. The 
eigenvalues provide a more precise picture of the way additional dimensions are explored, 
since they are a measure of the relative weight with which the dimensions of Hilbert space are 
explored. One sees a slow leaking of probability into additional dimensions. This leaking is 
due to the fact that, with an increasing number of perturbed steps, the probability increases 
for a state to have significant support outside the regular region, i.e., in the chaotic region. 
The part of the wave function that is in the chaotic region is subject to chaotic time evolution 
and therefore free to explore almost all dimensions of Hilbert space. 

VI. CONNECTION WITH STATISTICAL PHYSICS 

Consider a physical system, classical or quantum, with a known Hamiltonian. The state 
of the system at time t = to represents the observer's knowledge of the way the system 
was prepared. In classical physics, states are described mathematically by a probability 
density p(x) in phase space, while in quantum mechanics, states can be represented either 
by a Hilbert-space vector or, more generally, by a density operator p, depending on 
the preparation procedure. Observers with different knowledge assign different states to the 
system; the state is therefore not a property of the system alone, but refiects the observer's 
state of knowledge about the system. 

The entropy (in bits) of a system state, defined in the classical H = 

— J dTp{x.) log2[p(x)], where dT is the usual phase-space measure, and in the quantum case 
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as H = — Tr[plog2(p)], measures the information missing toward a complete specification of 
the system. The classical entropy is defined up to an arbitrary additive constant, reflecting 
the fact that an infinite amount of information would be needed to give the exact location 
of a single point in phase space. The quantum entropy vanishes for a pure state p = 
which is meaningful because no information beyond that contained in the wave function 
exists about a quantum system. As a consequence of Liouville's theorem, both classical and 
quantum entropy remain constant under Hamiltonian time evolution. 

To make the connection with thermodynamics, we assume that there is a heat reservoir 
at temperature Tq, to which all energy in the form of heat must eventually be transferred, 
possibly using intermediate steps such as storage at some lower temperature. In the presence 
of this fiducial heat reservoir, the free energy or maximum average extractable work for an 
equilibrium state is given by Fq = E—T^ks In 2 iJ, where E is the mean internal energy of the 
state. (More precisely, it is the difference between the Fq values of two states that determines 
how much work can be extracted in a transformation between the two states.) It is the main 
premise of this section that the maximum average extractable work, from now on called 
available work, is given by Fq = E — TQkBln2H for any state, even outside equilibrium. This 
means that each bit of missing information costs one TofcBln2 of available work. General 
arguments for this premise will be given elsewhere [0; here we discuss just one important 
case. If the non-equilibrium state was formed through reversible Hamiltonian time evolution 
starting from an initial equilibrium state, the amount of work Fq = E — TokBlii2 H can 
in principle be extracted on average if the system is made to evolve back into the initial 
equilibrium state using time reversal. This argument is equivalent to Loschmidt's famous 



Umkehreinwand Although time reversal appears to be impractical in most situations — 
a remarkable exception is spin echo — there are no known fundamental reasons for excluding 
it. 

Since entropy is a measure of the state of knowledge about the system, and since available 
work is determined by the entropy, the only way the available work can change (except for 
changes in the energy levels) is via a change in the state of knowledge. Hamiltonian time 
evolution of an isolated system does not lead to a change in the state of knowledge; entropy 
and available work remain unchanged. This is a consequence of Liouville's theorem and 
is true for regular as well as for chaotic systems: Hamiltonian time evolution in isolated 
chaotic systems does not lead to information loss The following paragraphs discuss 

the three ways in which information about the system can change: measurement, deliberate 
discarding of information, and interaction with an incompletely known environment. 

Available work can increase if an observation is made on the system. The accompanying 
decrease in entropy does not constitute a violation of the second law of thermodynamics, 
however, because the physical state of the observer changes in the course of the observation. 
Landauer ||2^, following seminal work by Szilard [^, has provided a simple and elegant 
quantitative description of the change in the state of the observer. If the observer wants to 
use additional information about the system to increase the available work, he must keep 
a physical record of the information. According to Landauer's principle, the erasure of a 
bit of information in the presence of a reservoir at temperature Tq is necessarily accompa- 
nied by dissipation of at least an amount To/cBln2 of energy. If this thermodynamic cost 
of erasing information is taken into account as a negative contribution to free energy, no 
observation can increase the total available work on the average P5|-p7[]. Information that 
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the observer possesses about the system therefore plays a role complementary to the en- 
tropy or missing information. Entropy or missing information about the system reduces the 
available work through the usual entropy term in the free energy; information the observer 
actually possesses must be stored physically, thus reducing further the available work due 
to the Landauer erasure cost. Total available work is determined by the sum of entropy and 
information [p5|-|27|J^. 

Available work can decrease if information about the system is lost. Information loss 
is equivalent to entropy increase. There are two main mechanisms leading to information 
loss: deliberate discarding of information and loss of information through interaction with 
an incompletely known environment. It must be emphasized that the well-known sensitivity 
to initial conditions in classical chaotic systems does not entail information loss, because 
statistical physics is concerned with the time evolution of probability distributions governed 
by Liouville's equation, not with trajectories of single phase-space points. 

Deliberate discarding of information was used by Jaynes to derive traditional 

thermodynamics. Jaynes showed how equilibrium thermodynamics follows effortlessly from 
Liouville's equation if only information about the values of the macroscopic variables defining 
a thermodynamic state is retained. In Jaynes's approach, irrelevant information is discarded 
by means of the principle of maximum entropy. Another way to discard information con- 
sidered irrelevant is coarse graining, where all details of a state below a certain scale are 
ignored. 

In contrast to these examples where information is discarded deliberately, an actual loss 
of information can occur in a system that, rather than being perfectly isolated, interacts 
with an incompletely known environment. The interaction with the environment leads to 
a perturbed time evolution of the system. Predictions for the system state are made by 
tracing out the environment, i.e., by averaging over the perturbations, which is generally 
accompanied by an entropy increase. 

Nothing forces one, however, to average over the perturbations. Alternatively one could, 
by making observations on the environment, gather enough information about the pertur- 
bations to keep track of the perturbed evolved system state to a certain accuracy, thereby 
reducing the entropy increase. In Section |V|, with the accuracy determined by the resolution 
angle 0, the minimum information needed to keep track of the system state to accuracy </> 
was denoted by AJ(0), and the resulting average entropy increase, the trade-off entropy, was 
denoted by AH{(f)). Averaging over the perturbations corresponds to an accuracy = 7r/2: 
AI{n/2) vanishes, and AH{7t/2) is the entropy due to averaging over the perturbation. 
Equation ( |4.12|) shows that the sum of information and trade-off entropy is never less than 
the entropy due to averaging, so that one can never gain in terms of total available work 
by gathering information about the perturbations. But at this stage, one has no reason 
to expect that one would do much worse by keeping track of the system state, so that in 
principle, the system entropy could be kept from increasing. 

For a system showing hypersensitivity to perturbation, however, there is a compelling 
reason not to keep track of a perturbed system state, but to average over the perturbations. 
The information needed to keep track of a perturbed state increases far more rapidly than 
the entropy due to averaging, which means that keeping track of the perturbed state would 
lead to an enormous reduction in available work due to the thermodynamic erasure cost. We 
have conjectured [0,0,0 that hypersensitivity to perturbation provides a quantitative link 
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between chaos and entropy increase in both classical and quantum open systems. Within 
the limits of our numerical method, the present paper establishes this link for the quantum 
kicked top. 
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APPENDIX A: VOLUME CONTAINED WITHIN A SPHERE IN HILBERT 

SPACE 

In this Appendix we compute the volume contained within a sphere of radius $ in D- 
dimensional Hilbert space. More precisely, we work in projective Hilbert space, i.e., the 
space of Hilbert-space rays or the space of normalized state vectors in which vectors that 
differ by a phase factor are equivalent. 

We begin by deriving the line element of the Riemannian metric that corresponds to 
the Hilbert-space angle ( |2.12[ ). Consider two neighboring normalized state vectors and 



+ \dilj). The infinitesimal angle ds between these vectors satisfies 

1 -rfs^ = cos^ds = |(^ I (1^) + \di)))'^ 

= l + 2Re((V^|#)) + |(V^|#)|2. (Al) 

The normalization of and + \dil)) implies that 

2Re((V'|#)) = -(#|#) , (A2) 

so the line element becomes 

ds^ = (#|#) - = {d^^\dxljj_) , (A3) 

where \dip±) = {dip) — lip) {ipldip) is the projection of orthogonal to \ip). The metric ( |A3[) , 



called the Fubini-Study metric |^, is the natural metric on projective Hilbert space. Notice 
that the line element is invariant under phase changes of {ip) and {ip) + {dip). 

Consider now a sphere of radius $ < tc/2 in projective Hilbert space; let the center of 
the sphere be denoted by \ipo). Any normalized vector \ip) can be written as 

\^) = e"-^ cos (pltpo) + sin (f)\r]) , (A4) 

where |?7) is a normalized vector in the subspace orthogonal to \ipo), and the polar angle 
(f) = cos~^ {\{ip\ipo)\) satisfies < < 7r/2. The region contained within our sphere of radius 
$, which we denote by Vd{^), consists of all vectors such that < $. The phase freedom in 
\tp) can be removed by choosing the phase 6 = 0. That having been done, \ri) ranges over all 
normalized vectors in the subspace orthogonal to iV'o); iii particular, two normalized vectors 
|?7), differing only by a phase factor, are not equivalent. 
We can now write 
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Idip) = — sin (p d(f)\ipo) + cos (f> d(f)\r]) + sin (p\dri) , (A5) 

where \dri) is the infinitesimal change in \ri) (notice that (ipoldf]) = 0). Normahzation of |?7) 
and \ri) + \dri) imphes, just as in Eq. (|A2|) , that 

2Re{{ri\dri)) = -{dr]\dr]) . (A6) 

Retaining only second-order terms in the infinitesimal changes, we can put the line ele- 
ment (|A3|) in the form 



ds^ = dcj)"^ + sm^ (t)d^'^ , (A7) 

where 

d-i'^ = {dT]\dT]) - siYi^ <p\{T]\dr])\'^ (A8) 

defines a Riemannian metric on the the space of normalized vectors in the l)-dimensional 
subspace orthogonal to lipo). This space is a {2D — 3)-dimensional sphere of unit radius, 
denoted S2D-3, although the line element is, as we show below, different from the 
standard geometry on a {2D — 3)-dimensional unit sphere. Notice that the Fubini-Study 
metric ( [A7| ) scales all lengths on S2D-3 by a factor of sin 0; this scaling is analogous to the 
way that polar angle on an ordinary 2-sphere scales the size of circles (1-spheres) of latitude. 

Consider now any orthonormal basis \r]j), j = 1, D — 1, in the subspace orthogonal 
to \ipo). We can introduce coordinates on S2D-3 by expanding |?7) as 

D-l 

\V) = T.i^J+wM) ■ (A9) 
i=i 

Normalization of \ri) implies the constraint 

D-l 

which defines the {2D — 3)-dimensional unit sphere. The first term in the metric (|A8|), 

D-l 

{dri\dri) = J2 dx] + dy] , (All) 

i=i 

is the fiat Euclidean metric in 2 (D — 1) dimensions; it induces the standard metric on 5*20-3 • 
The second term in the metric (|A8|) modifies the standard geometry on S2D-3- Noting that 



'D-l \ D-l D-l 



{rj\dri) = ]:dl J2 ^'j + ] + ^ J2 ^j^yj - yjdxj = i ^ Xjdyj - yjdxj , (A12) 
we can put the second term in (|A8|) in the form 



'D-l ^ 2 



sin^ (f) \ {'r]\dr])\'^ = sin^ ^ Xjdyj — yjdxj . (A13) 
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Perhaps the easiest way to see how the second term affects the geometry on S2D-3 is to 
make a judicious choice of coordinates. Given an arbitrary vector \ri), we can always choose 
the orthonormal basis so that |?7i) = \ri), which means that \ri) is assigned coordinates Xi = 1, 
Hi = 0, and Xj = yj = 0, j = 2, . . . , D — 1. In these special coordinates the metric ( |A8D , 
evaluated at \ri), takes the form 

D-l 

d'j'^ = cos^ (t)dyl+ dx] + , (A 14) 

where we have used the constraint ( |A1CI|) to write 



D-l 



dxi = — I —yidyi — ^ Xjdxj + yjdyj ] =0 aX \r]). (A15) 

i=2 



Xi 



In these special coordinates the standard geometry on 820-3 has the same form at \ri), 
except that the cos0 is replaced by 1. The effect of the second term in the metric ( |A8D is 
thus to shorten lengths (relative to the standard geometry) along one direction on S2D-3 by 
a factor cos0; the direction of shortened length corresponds, in the special coordinates, to 
the yi direction or, in coordinate-free language, to an infinitesimal change in the phase of 

In the special coordinates the volume element on S2D-3 at |?7), defined by the line element 
d'j'^, is given by 

cos dyidx2 ■ ■ ■ dxD-idy2 ■ ■ ■ dya-i = cos dS2D-3 , (A16) 

where dS2D-3 is the volume element defined by the standard geometry on 820-3- Writing 
the volume element ( |A16D in terms of dS2D-3 frees it from dependence on the special coor- 
dinates. Referring to the Fubini-Study metric ([^), we can now write the volume element 
on projective Hilbert space as 

dVo = (sin 0)2^-=^ cos c/0 rf^sD-s • (A17) 

The 2D — 3 factors of sin0 come from scaling all 2D — 3 dimensions of S2D-3, as required 
by the Fubini-Study metric ([A7| ). 

We are now prepared to compute the volume Vz)($) contained within a Hilbert-space 
sphere of radius $: 



V^($) = / dVD= I d0 (sin 0)2^-3 cos / dS2D-3 
Jvnm Jo J 

(sin $)2(^-i) = (sin $)2(^-^) . (A18) 



<S2D-3 

2(/^-l) 



Here S2D-3 is the volume of the unit sphere S2D-3, calculated using the standard geometry 
(Be careful: this is the "area" of 6*2^-3, not the volume interior to it), and 

Vn = Vn{n/2) = (A19) 
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is the total volume of projective Hilbert space. 
The volume of an n-dimensional unit sphere S'„, 



27r("+i)/ 

r 



= ,^^1, , (A20) 



2 

follows from a standard trick involving Gaussian integrals: 

oo „ / /-oo „\ n+1 



drr'^Sne-'" = { due^""") = . (A21) 



This gives us S2d~3 = 27r ^/ — 2)!, which allows us to write the total volume of projective 



Hilbert space ^1 



as 

= ■ 

We can obtain immediately two other results: (i) the volume contained between two 
Hilbert-space spheres of radius (j) and (j) + dcj) is 

(iV^(0) = 52D-3(sin0)2^^^cos0# = Vd2{D - l){sm (f))^^-^ cos (f) d(f) ; (A23) 

(ii) the probability that two vectors selected at random are separated by a Hilbert-space 
angle between and (f) + d(f) is 

g{(j))d(j)= '^\''}'^^ = 2(D -1) (sin 0)^^-=^ cos 06/0. (A24) 



APPENDIX B: ENTROPIES OF DISTRIBUTIONS OF VECTORS WITHIN A 

HILBERT-SPACE SPHERE 

In this Appendix we compute entropies of density operators constructed from vectors 
that lie within a sphere of radius $ in D-dimensional Hilbert space. We first compute the 
entropy Ho{^) of a density operator p that corresponds to a uniform distribution of Hilbert- 
space vectors in the region Vo{^) contained within a sphere of radius Formally, p is given 
by 

P= f , ,^^I^)(V^I, (Bl) 



where Vz)($) = (sin ^Y^^^^^Vd [Eq- (|A18|) ] is the volume contained within a sphere of radius 
$ in D-dimensional Hilbert space, and dVo = {sin cos (j) dcp dS2D-3 [Eq- ( |A17D ] is the 

volume element on projective Hilbert space. Let the center of the sphere be denoted by \ipo). 
Symmetry about l-^o) entails that p have the form 
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P = Ao + (i - |^o)(^o|) . (B2) 



The eigenvalues of p are 



dVo 1/ / I / \|2 [ ^Vd 2 



Ao = (^olpl^o) = J^^^^^ K^o|^)| = J^^^^^ ^ cos 



' ''O SlIKB COS < 



(sin $)2(^-i)Vd ^ 



and 



- — - sin' $ > - (B3) 
D-D ^ ' 



D-1 D - D 
The entropy of p is thus 



i/,,($) = -5:A,.log2A 



fe=0 



-Ao log2 Ao - (1 - Ao) log2(l - Ao) + (1 - Ao) hg.iD - 1) (B5) 
sm $ - 1 logo 1 sni' $ 



D / D 

D-l o ^ , /sin2$\ 
-_sin2$log2f-^j . (B6) 

For comparison, the spread of p [Eq. ([4.7D] has its maximum value, regardless of the value 
of $, 

H2A^) = -''p^Y^Jog,(^-^^=\og,{D-l), (B7) 

reflecting the fact that the uniform distribution of vectors within a sphere of radius $ explores 
all D — 1 dimensions in the subspace orthogonal to {ipo)- This result regarding the spread does 
not depend on the particular value of Ao: any density operator of the symmetric form (B2) 
has the maximum spread, log2(-D — 1), provided only that Ao is the largest eigenvalue, i.e., 
Ao > l/D [cf. Eqs. (pD and (lESl)]. 

It is interesting to compare the entropy of a uniform distribution of vectors in Vd(<I') 
with the maximum entropy that can be attained by distributing vectors in Vr>{^). To do 
this, consider a density operator 

N 
i=l 

constructed from vectors \ipi) that lie in Vd^^), i.e., | (V'olV'i) I = cos0i > cos$, i = 1, . . . , N. 

A unitary transformation in the subspace orthogonal to |?/^o) leaves the angles 0j and 
the entropy unchanged. Mixing the density operators that result from all such unitary 
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transformations gives a new density operator which is symmetric about \ipo), but which, 



because of the concavity of the entropy [p!9iP0i , has an entropy that is not smaller than the 



entropy of p. Thus, in seeking the maximum entropy, we can restrict attention to density 
operators that are symmetric about {i/jq) and, hence, have the form ( P2|) and have entropy 
given by Eq. ( |B5D . This entropy, which is a function of the eigenvalue Aq, has a maximum 
value of log2 D at Xq = 1/D and decreases monotonically away from this maximum in either 
direction. 

The eigenvalue Aq is bounded below by 

N 

Ao = (V^oIpIV^o) = Y.Pi cos^ ^ cos^ ^ • (B9) 

i=l 

Hence, an upper bound on the entropy follows from choosing Aq = cos^ $ when cos^ $ > 1/D 
and choosing Aq = l/D when cos^ $ < The upper bound is given by 

.max 

($) 

_ J - cos^ $ log2 cos^ $ - sin^ $ logs sin^ ^ + sin^ $ log2(-D - 1) , cos^ $ > l/D , , . 
~\\og^D, cos2$<l/D. ^ 

This function is plotted in Fig. ^ for D = 512. Notice that the entropy ( |B6| ) of a uniform 
distribution of vectors within a sphere of radius $ approaches the upper bound ( |BlUp as 
D ^ oo. 

That the bound is actually the maximum, as implied by the notation, is demonstrated by 
finding a density operator, constructed from vectors in V/:)($), which achieves the entropy 
upper bound. To that end, consider the 2{D — 1) vectors 

= cos<l>\iJo) + sm<!>\r]j) , j = 1,...,D -1 , 
|7/;;) = cos$|V'o)-sin$|r7,), j = l,...,D-l, (Bll) 

where the vectors \7]j) make up an orthonormal basis in the subspace orthogonal to {ipo). 
The vectors (|B11| ) all lie on the sphere of radius $ — as far from \tpQ) as is allowed. Construct 
the density operator 



D-l 



= (cos^ $ + posin^ $)|V'o)(^o| + -^^^ sin^ $ (1 - |^o)(^o|) , (B12) 

which has the symmetric form of Eq. ( p^) , with Aq = cos^ $+Po sin^ $. One would obtain the 
density operator ([B12|) by letting p— Pol'^o) (V'ol be constructed from any set of vectors that lie 



on the sphere of radius $ and are symmetrically distributed about {ipo). To achieve the upper 
bound ( [B10| ), one chooses po = when cos^ $ > l/D and chooses po = (1/-D — cos^ $)/ sin^ $ 
when cos^ $ < l/D. 
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FIGURES 



FIG. 1. Distribution of 2" vectors randomly chosen in i^-dimensional Hilbert space. Each 
diagram shows, as a function of the angle (j), the distribution g{(j)) of Hilbert-space angles (unnor- 
malized, in arbitrary units), the average information A/(0) to specify a vector at the resolution 
given by (j) (in bits), the trade-off entropy AH{(j)) (in bits), and the average spread AH2{(f)) (in 
bits). For a precise definition of these quantities, see Section (a) n = 12, D = 512. (b) n = 10, 
D = 62. 



FIG. 2. Results characterizing the distribution of Hilbert-space vectors for the perturbed kicked 
top with J = 511.5 and k = 3 m the presence of the twist perturbation (3.1) with e = 0.03. The 
same quantities as in Fig. |l] are shown, (a) Chaotic case, i.e., initial coherent state \0,ip) centered 
in the chaotic region with 9 and given by Eq. (p. 7]). Distribution of all 2^ vectors generated after 
8 perturbed steps, (b) Regular case, i.e., initial coherent state centered at the elliptic fixed point 
given by Eq. ( ^.51) . All 2^ vectors generated after 8 perturbed steps, (c) Chaotic case. 12 steps, all 
2^^ vectors, (d) Regular case. 12 steps, all 2^"^ vectors. 



FIG. 3. Eigenvalues of the density operators formed by averaging over random vectors 
(squares), vectors generated by the perturbed kicked top in the chaotic case (diamonds), and 
vectors generated by the perturbed kicked top in the regular case (crosses). The density operator 
in the random case was generated by averaging over 1024 vectors randomly chosen in 62-dimen- 
sional Hilbert space. The density operators in the chaotic and regular cases were formed from the 
2^^ vectors in Figs. ||(c) and (d), respectively. The main diagram shows, for all three cases, the 
largest 62 eigenvalues greater than 10"^*^. The inset shows all eigenvalues greater than 10~^^ for 
the chaotic case (solid line) and the regular case (dashed line). 



FIG. 4. As Fig. ^, but using the twist perturbation (|3.1| ) with e = 0.003. (a) Chaotic case. 12 
steps, all 2^^ vectors, (b) Regular case. 12 steps, all 2^"^ vectors, (c) Chaotic case. 2^^ vectors 
randomly chosen after 200 perturbed steps, (d) Regular case. 2^*^ vectors randomly chosen after 
200 perturbed steps. 



FIG. 5. Upper bounds for the average information, A/£)_jnax! and for the trade-off entropy, 
-f^D,max, for 2^'^'^ random vectors in 512-dimensional Hilbert space, as given by Eqs. ( ^.21) and (|B10D . 
The curve for Al/j max is an excellent approximation to the average information A/((/>) except for 
a small region around (p^. The trade-off entropy AH{(f)) is well approximated by the curve for 
HD,ma.x{<t>) for angles above 0;,, but goes to zero rapidly for angles below (f)},. 



FIG. 6. As Fig. |, but using the turn perturbation ( ^ ) with e = 0.003. (a) Chaotic case. 12 
steps, all 2^^ vectors, (b) Regular case. 12 steps, all 2^^ vectors, (c) Chaotic case. 2^^^ vectors 
randomly chosen after 30 perturbed steps, (d) Regular case. 2^*^ vectors randomly chosen after 30 
perturbed steps. 
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FIG. 7. As Fig. ^, regular case, (a) 2^^ vectors randomly chosen after 100 perturbed steps, 
(b) 2^^^ vectors randomly chosen after 200 perturbed steps, (c) The 20 largest eigenvalues of the 
density operators obtained by averaging over 2^*^ vectors randomly chosen after 30, 100, and 200 
perturbed steps. 
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